Correlated defects, metal-insulator transition, and magnetic order 
in ferromagnetic semiconductors 
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The effect of disorder on transport and magnetization in ferromagnetic III-V semiconductors, 
in particular (Ga,Mn)As, is studied theoretically. We show that Coulomb-induced correlations of 
the defect positions are crucial for the transport and magnetic properties of these highly compen- 
sated materials. We employ Monte Carlo simulations to obtain the correlated defect distributions. 
Exact diagonalization gives reasonable results for the spectrum of valence-band holes and the metal- 
insulator transition only for correlated disorder. Finally, we show that the mean-field magnetization 
also depends crucially on defect correlations. 

PACS numbers: 75.50.Pp, 71.30.+h, 72.20.Ee 
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Introduction. — Recently, there has been substantial in- 
terest in diluted ferromagnetic III-V semiconductors due 
to observations of relatively high Curie temperatures 
iQ, H, ^. This makes these materials promising for ap- 
plications as well as interesting from the physics point of 
view. They could allow the incorporation of ferromag- 
netic elements into semiconductor devices, and thus the 
integration of processing and magnetic storage on a single 
chip. To be specific, we consider here manganese-doped 
GaAs. The properties of this material rely on the dual 
role played by the Mn impurities: They carry a local spin 
due to their half-filled d-shell and dope the system with 
holes, which mediate a ferromagnetic indirect exchange 
interaction between the spins. Furthermore, the mate- 
rials are highly compensated presumably due to antisite 
defects (As substituted for Ga) [Q, which drives them 
towards a metal-insulator transition (MIT). It is clearly 
important to understand the interplay between transport 
properties, magnetic ordering, and the defect configura- 
tion. In this letter we show that correlated defects are 
required for a description consistent with experiments. 

There is a growing body of theoretical work on 
(Ga,Mn)As §. Most approaches either start from the 
heavily doped regime P, 0, ^, ||, |lO[ |ll|, |l^ or from 
the weak-doping limit pj^ . In the case of heavy dop- 
ing ^ H, H, |l^, pr], [12[ one considers valence-band 
holes moving in the disorder potential of the defects. For 
large hole concentration the Fermi energy Ep lies deep 
in the valence band compared to the Coulomb-potential 
fiuctuations and the latter are neglected. On the other 
hand, for light doping the local impurity states over- 
lap only weakly, forming an impurity band. This model 
is problematic since for the typical doping range the im- 
purity band would be much broader than the energy gap 
[Q . Both approaches have in common that the effect of 
antisites on the disorder and correlations of defects are 
neglected. 

The outline of the remainder of this paper is as follows. 
First, we discuss the configuration of Mn impurities and 
antisite defects with the help of Monte Carlo simulations 



and obtain the disorder potential and its spatial corre- 
lations. Then, we derive the spectrum and localization 
properties of holes in this potential and discuss the MIT. 
Finally, we calculate the polarizations of Mn and hole 
spins within a selfconsistent mean-field theory. 

Defect configurations. — In the theory of doped semi- 
conductors one often assumes the dopants to be randomly 
distributed ||^, ^ . However, the present system is highly 
compensated. This leads to defects of either charge be- 
ing present (relative to Ga, Mn impurities have charge 
— e while antisites carry +2e) and to a low density of 
carriers, which only weakly screen the Coulomb interac- 
tion between the charged defects In the resulting 
Coulomb plasma the defect positions are expected to be 
highly correlated. We study here the limiting case where 
the defects come into thermal equilibrium during growth 
I p^ . This is motivated by experiments Q suggesting 
that defects diffuse rapidly at 250°C, which is a typical 
growth and annealing temperature for (Ga,Mn)As grown 
by molecular beam epitaxy. The resulting configuration 
is assumed to be quenched at low temperatures. 

To find typical defect configurations close to equi- 
librium we perform Monte Carlo simulations for the 
classical Coulomb system of Mn impurities and anti- 
sites on the Ga sublattice with the Hamiltonian H = 
1/2 ^j^- (7igj/(ery) e~'''j/''°", where qt are the defect 
charges, is their separation, and e = 11. The screen- 
ing length Tscr is obtained from nonlinear screening the- 
ory p5[ |. Tscr is much larger than the nearest-neighbor 
separation of Ga sites for realistic parameters so that it 
hardly affects the small-scale defect correlations relevant 
here. We employ the Metropolis algorithm at 250°C for 
systems of 20 x 20 x 20 conventional cubic unit cells with 
periodic boundary conditions, unless stated otherwise. 

A typical result for Mn concentration x = 0.05 and 
p — 0.3 holes per Mn, as suggested by experiments [Q, 
is shown in Fig. |^. The concentration of antisites is ob- 
tained from X and p under the assumption of charge neu- 
trality. Most of the defects have formed clusters. This 
leads to screening of the disorder potential by the de- 
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FIG. 1: Typical equilibrium configuration at 250° C for 5% 
Mn doping and 0.3 holes per Mn. Only the defects are shown, 
large (small) circles denote antisite (Mn) defects. Defects at 
nearest neighbor Ga sites are connected by a bond. For clarity 
we show only 10x10x10 conventional unit cells. 
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FIG. 2: Potential correlation function D{r) for 5% Mn and 
0.3 holes per Mn, both for uncorrelated defects and for equi- 
librium at the growth temperature. The nearest-neighbor Ga 
site separation of 3.99 A is also indicated. 



fects, as can be seen from the potential correlation func- 
tion D{r) = {V{r) V^(r'))|,-,'|=, - {V)^ plotted in Fig. | 
for X = 0.05 and p = 0.3. Note that AV = y/D{0) 
is the width of the distribution of V'(r). For correlated 
defects, D{r) and thus AT^ are strongly reduced and the 
Coulomb potential is screened on the much shorter length 
scale Tion, which is of the order of the nearest-neighbor 
Ga site separation (This ionic screening is only 

possible due to the presence of defects of either charge. 
Thus it does not take place, e.g., in ferromagnetic II- VI 
semiconductors since Mn is isovalent in these materials.) 
Importantly, however, AT^ w 0.284 eV is still not small 
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FIG. 3: IPR as a function of energy at the valence-band edge 
for 5% Mn and 0.3 holes per Mn. The thin (heavy) solid line 
denotes empty states for uncorrelated (correlated) defects, the 
dotted lines denote occupied states. 



compared to the Fermi energy \Ep \ « 0.329 eV. 

Recent high resolution x-ray diffraction [|l8j as well as 
resistivity and magnetization measurements support 
the formation of small clusters. The more direct determi- 
nation of the local structure around Mn dopants by ex- 
tended x-ray-absorption fine structure ]l9| does not give 
direct information on the abundance of specific elements 
around the Mn. However, the observed loss of order 
around Mn is certainly consistent with cluster formation. 
Finally, we show in the following that experimental re- 
sults for the band gap, transport, and the magnetization 

^ strongly support the clustering. 

Hole states. — To find the spectrum and localization 
properties of valence-band holes moving in the disor- 
der potential V{r), we start from the Hamiltonian H = 
-E.(?iV2m*) V| + F(r,), where we use the envelope 
function and parabolic-band approximations The 
calculations are done for spin-less holes, which is justified 
since the additional disorder introduced by the exchange 
interaction with the Mn spins is much smaller than AV. 

The Hamiltonian is solved by exact diagonalization in 
a plane- wave basis of IS'^ = 2197 states giving the 
energy spectrum and normalized eigenfunctions. The 
Fermi energy Ep is obtained from the hole concentra- 
tion px. We also obtain the inverse participation ra- 
tio IPR(n) = 1/Er \'Pnir)\^ of the states V™(r). The 
IPR allows to estimate the position of the mobility edge, 
since it is of the order of the volume for extended states 
but is essentially independent of system size for local- 
ized states. Figure || shows a comparison of the IPR 
as a function of energy for uncorrelated defects and for 
the equilibrium configuration obtained above, for 5% Mn 
and 0.3 holes per Mn. The valence-band edge is strongly 
smeared out by uncorrelated defects so that the energy 
gap is completely filled, which is not observed experimen- 
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FIG. 4: IPR as a function of energy at the valence-band edge 
for various Mn dopings and hole concentrations p per Mn 
taken from experiment 111, ESl. 



tally, whereas correlated defects only lead to a small tail 
in the gap. We have also studied the dependence of the 
IPR on system size over the limited feasible range in 
L, finding that IPR(n) cx for the flat part of the curve, 
signifying extended states, whereas the IPR is indepen- 
dent of L for states sufficiently far in the tail, which are 
thus localized. Hence, the mobility edge lies on the slope 
of the curve. The states are much more localized for un- 
correlated defects, due to the larger disorder potential, 
whereas for correlated defects the states close to Ep are 
clearly extended. Thus the spectrum and the localization 
properties agree with experiments only for correlated 
defects. 

In Fig. ^ we show the IPR as a function of energy for 
correlated defects at various Mn dopings and hole concen- 
trations. The hole concentration for given Mn doping was 
chosen in accordance with experiments [Q, |23j . The spec- 
trum and the IPR change little with Mn doping, since 
the potential screened by the defects on the length scale 
Tion is much less affected by a change of concentration 
than the bare Coulomb potential. 

The dominant effect comes from the hole concentra- 
tion. Our result for the scaling of the IPR suggests that 
for 5% and 3% Mn the states at Ep are extended and we 
expect metallic behavior. For 1% Mn Ep is close to the 
mobility edge, while for 0.5% Mn the states at Ep are 
localized. From this we estimate the MIT to take place 
at a Mn concentration of the order of 1%. This result is 
in reasonable agreement with experiments Q|. For com- 
parison, uncorrelated defects would lead to an MIT at 



about 5% Mn, see Fig. |. 

Magnetization. — Finally, we discuss the spontaneous 
magnetization and its dependence of the defect con- 
figuration. The exchange interaction between hole 
and Mn spins is described by the Hamiltonian H = 

where V{r) is the disorder potential obtained above, 
-'^pd ~ -45eVA3 is the exchange integral [|[ |4[ H), 
Si and S; are the hole and Mn spin operators, respec- 
tively, and R; are the Mn impurity positions. In second 
quantized form, the Hamiltonian reads 

--^Pd "^'"^ V4(Ri) ■ Si 1pn'(R-l) Cn'a', 

where e„ are the eigenenergies obtained in the absence of 
exchange as discussed above, ipnii^) are the correspond- 
ing eigenfunctions, and r is the vector of Pauli matri- 
ces. The exchange interaction is decoupled at the mean- 
field level, introducing the averaged Mn spin polariza- 
tions M; = (S;) and hole spin polarizations at the Mn 
sites, m; = (i;„„,„v'cL'/';(R/)(w/2)?/'„'(R/)c„,<,-). 
We do not perform a spatial average, i.e., the disorder 
is retained. The hole and Mn sectors can now be diag- 
onalized separately, M/ and m/ are calculated, and the 
procedure is iterated. This selfconsistent mean-field the- 
ory is similar to the one employed by Berciu and Bhatt 
p3[ . The main difference is that we start from realistic 
hole states for the disorder potential. 

Figure ^(a) shows the Mn and hole spin polarizations 
for X — 0.05 and p — 0.3. Also shown are the polariza- 
tions obtained neglecting disorder by using plane waves 
for the hole states and the virtual crystal approximation 
(VGA) for the Mn spins ||, |, |2l). We see that disorder 
strongly enhances Tc jl^. This increase is easily under- 
stood: Tc is determined by the indirect exchange inter- 
action of nearest-neighbor Mn pairs, which is inversely 
proportional to the separation 0, ]|]. Without disor- 
der, the typical distance is r^ca ^ "mIi ' where tt-mii is the 
density of Mn impurites. For x = 0.05, Tvca 9.67 A. On 
the other hand, with disorder the nearest-neighbor sep- 
aration on the Ga sublattice is rdis = 3.99 A. Therefore, 
the indirect exchange, and thus Tc, is larger by a factor 
of order rvca/^dis ~ 2.4 in the disordered case. 

For correlated defects the polarization curve looks 
more "mean-field-like," in qualitative agreement with ex- 
periments ||l|, |l^, 1^. This can be explained by the ex- 
tended nature of the relevant hole states for correlated 
defects, which leads to a long-range indirect exchange in- 
teraction. Thus the effective field seen by a Mn spin is 
averaged over many other Mn spins and the spatial fluc- 
tuations of this effective field are small. On the other 
hand, for uncorrelated defects the holes are partly local- 
ized, the indirect exchange is of shorter range, and the 
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FIG. 5: (a) Magnetization as a function of temperature for 5% 
Mn and 0.3 holes per Mn, both for uncorrelated and correlated 
defects. The results denoted by "VGA" neglect disorder, (b) 
Magnetization for various Mn and hole concentrations (with 
correlated defects). 



fluctuations are larger. We believe that this leads to the 
smeared-out magnetization curve. It is clearly important 
to study the effect of fluctuations on the magnetization. 
We expect fluctuations to reduce Tc in particular for the 
uncorrelated case. Figure ||(b) shows that is reduced 
for smaller Mn and hole concentrations and that the typ- 
ical temperature scale is approximately proportional to 
X, in accordance with experiments Note also that in 
all cases the holes are only partially polarized. 

To conclude, we have obtained typical defect config- 
urations in (Ga,Mn)As and studied the effect of the re- 
sulting disorder potential on the valence-band holes. The 
strong Coulomb interactions of charged defects leads to 
the formation of defect clusters during growth. The pre- 
viously neglected antisite defects are crucial for this ef- 
fect. We have shown that such a correlated defect distri- 
bution is required to understand the hole spectrum, the 
metal-insulator transition, and the shape of the magne- 
tization curves. Our results should also help to clarify 
why the properties of ferromagnetic semiconductors de- 
pend so strongly on details of the growth process. — We 
profited from discussions with P. J. Jensen, J. Konig, J. 
Schliemann, and in particular M. E. Raikh, who empha- 
sized to us the importance of compensation. 

Note added. — Recently, Yang and MacDonald have 
studied the effect of the hole-hole interaction, assuming 
a random distribution of Mn dopants and antisites. 
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